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SUMJ4ARI 

The  purpose  of  this  analysis  Is  to  investigate  the  use  of  an 
automatic  control  to  damp  undesirable  lateral  oscillations.  The  type  of 
control  under  consideration  is  a  yaw  damper  servomechanism  with  rudder 
control  applied  in  proportion  to  yawing  velocity.  The  effect  of  angle 
of  attack  and  gyro  axis  angle  on  components  of  roll  and  yaw  picked  up  is 
developed,  and  the  optimum  angle  of  tilt  of  the  gyro  axis  is  selected. 
The  problems  of  maintaining  equal  effectiveness  of  the  control  through  the 
speed  range,  eliminating  yaw  compensation  during  steady  turning  or  rolling 
maneuvers,  critical  time  lag,  and  prevention  of  feedback  to  the  pilot's 
controls  are  considered. 

The  airplane  used  for  the  analysis  is  an  operational  Jet  fighter 
of  hig^  performance  which  has  been  reported  to  have -«  "snaking,"  or  "Dutch 
Roll"  problem.  The  flight  conditions  investigated  are  cruising  at  35,000 
feet,  50^  normal  rated  power  at  10,000  feet,  and  power  approach  at  sea 
level.  Calculations  in  these  configurations  predict  that  the  aircraft  has 
an  unsatisfactory  "Dutch  Rolling"  tendency,  probably  in  all  fligjit  condi- 
tions. Installation  of  the  proposed  damper  system  will  enable  the  airplane 
to  satisfactorily  meet  the  lateral-directional  oscillation  requirements. 
In  addition,  the  stabilizer  will  cause  the  aircraft  to  be  spirally  stable, 
and  will  have  the  effect  of  gearing  the  rudder  to  the  ailerons  for  improved 
entry  into  rolling  or  turning  maneuvers. 
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SYMBOLS  AND  COEFFICIENTS 

a  angle  of  attack  meaaiired  from  the  flight  path  to  the  airplane 

reference  axis,  degrees 

6  angle  between  the  reference  axis  and  the  horizontal  axis,  positive 

when  the  reference  axis  is  above  the  horizontal  axis,  degrees 

"iJ  angle  of  the  flight  path  to  the  horizontal,  positive  in  a  climb, 

degrees 

y\  angle  of  attack  of  the  principal  longitudinal  axis  of  the  airplane, 

positive  when  the  principal  axis  is  above  the  flight  path  at  the 
nose,  degrees 

c  angle  between  the  reference  axis  and  the  principal  axis,  positive 

when  the  reference  axis  is  above  the  principal  axis,  degrees 

l)  angle  between  the  gyro  spin  axis  and  the  reference  axiSf  positive 

whysn  the  reference  axis  is  above  the  gyro  axis,  degrees 

k  angle  between  the  flight  path  and  the  gyro  axis,  positive  when  the 

gyro  axis  is  above  the  fli^t  path,  degrees 

V  angle  of  axirauth,  radians 
0  angle  of  bank,  radians 

P  angle  of  sideslip,  radians  (v/V) 

V  sideslip  velocity  along  the  Y  axis,   feet  per  second 

V  true  airspeed,   feet  per  second 

p  mass  density  of  air,  slugs  per  cubic  feet 

<r  density  ratio,  p/p 

q  dynamic  pressure,  pounds  per  square  foot  (l/2  pV2) 

W  wei^t  of  airplane  in  pounds 

m  mass  of  airplane  in  slugs,  W/g 

g  acceleration  of  gravity,  feet  per  second  per  second 

M  relative  density  factor,  n/pSb 

m  factor  for  dimensionallzation,  2b  M-A  " 

Cj^  trim  lift  coefficient,  W  cos)f/q3 


PRINCETON  UNIVERSITY  PAGE 

AERONAUTICAL  ENGINEERING  LABORATORY  report     tiz 


M  niach  number 

b  wing  span  in  feet 

Jl.         distance  from  the  center  of  gravity  of  the  airplane  to  the  center 
of  pressure  of  the  fiJi,   feet 

S  wing  area  in  square  feet 

S  control  surface  area  for  yaw  damper  in  square  feet 

S  .       vertical  tail  area  in  square  feet 

S  rudder  area  in  square  feet 

h  distance  from  the  longitudinal  body  axis  to  the  center  of  pressure 

of  the  vertical  control  surface 

D  the  operator,  d/dt 

P  period  of  oscillation  in  seconds 

T,  /„     time  for  amplitude  of  oscillation  to  damp  to  one-half  of  its 
'        original  amplitude  in  seconds 

C,  /„     cycles  for  the  amplitude  of  oscillation  to  damp  to  one-half  of  its 
'        original  amplitude 

K  control  gearing,  or  gain,    |    \^; 

K  amplitude  ratio  of  damper,     l/K  or      I    /c        n 

c  *      »       /  1    '  o  ri   Damper 

K.         airplane  amplitude,      |  ^/c     \ 

I  moment  of  inertia  about  the  principal  longitudinal  axis  in  slug 

o       feet  squared 

I  monient  of  inertia  about  the  principal  normal  axis  in  slug  feet  squared 

o 

I'  moment  of  inertia  coefficient  about  the  principal  longitudinal  axis, 

^         I^  /qSb 
o 

I  moment  of  inertia  coefficient  about  the  fli^^t-path  axis, 

(I'   cos^'n    -f  I'   sln^n   ) 

I  moment  of  inertia  coefficient  about  axis  normal  to  fli^t  path, 

(!•   cos  Ti     +  I'   sin  'n     ) 

Z  X  ' 

I'    moment  of  inertia  coefficient  about  the  principal  normal  axis.  I  /qSb 
z  z 

o 
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I  product  of  inertia  coefficient  with  respect  to  flight  path  axis  and 

^  axis  normal  to  flight  path,       -(I'   -   iMsin-^  cos-^) 

C^  rolling  moment  coefficient  (£2UiS^^on«ni) 

C_  yawing  moment  coefficient  (^ ^rr ) 

n  '    "  qbD 

C,  side  force  coefficient  (^^^.era^  force ^ 

r,H^  angular  yawing  velocity  in  radians  per  second 

p,  0  angular  rolling  velocity  in  radians  per  second 

dCi  dC 
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dp 
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The  stfimdard  NACA  conventions  in  regard  to  directions  of  axes»  angles, 
control  deflections,  forces,  moments,  and  linear  and  angular  motions  are  used. 
Where  dotted  quantities  are  written,  they  indicate  time  derivatives, ie,f  «  ^  » 
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INTRODUCTION 

The  purpose  of  this  report  is  to  conduct  a  theoretical  investigation 
into  the  use  of  a  yaw  damper  servomechanism  to  in5)rove  the  lateral  dynamics 
of  an  operational  high  speed  jet  filter.  Calculations,  flight  tests,  and 
pilots  reports  indicate  that  this  aircraft  does  not  satisfactorily  meet  the 
requirements  of  NAVAER  SR-li9B,  "Specifications  For  Flying  Qualities  of 
Piloted  Airplanes,"  in  regard  to  the  damping  of  the  lateral-directional  os- 
cillation, controls  free  or  fixed. 

High  performance  requirements  of  present  day  aircraft  and  those  of 
the  immediate  future  are  leading  to  designs  in  which  lateral  dynamic  flight 
stability  is  becoming  a  serious  problem.  "Dutch  roll"  at  low  airspeeds, 
particularly  on  landing  approaches,  is  annoying,  fatiguing,  and  dangerous 
to  the  pilot.   It  makes  instrument  approaches  particularly  difficult.  The 
present  day  trend  toward  sweepback  and  low  aspect  ratio  causes  high  effective 
dihedral  at  low  speeds  and  also  reduces  the  damping  in  roll  and  increases  the 
adverse  yaw  due  to  roll.  These  factors  combine  to  contribute  to  unsatis- 
factory lateral  dynamics.  "Snaking,"  the  control  free  lateral  dynamic  oscil- 
lation, compromises  the  effectiveness  of  the  aircraft  in  carrying  out  its 
mission  as  a  gunnery  or  bombing  platform.  Satisfactory  solutions  to  the 
problems  of  "Dutch  roll,"  and  "snaking"  will  not  only  improve  the  flying  qual- 
ities, but  will  give  the  designers  more  freedom  in  still  hi^er  performfiuice 
designs. 

Considerable  work  ia  being  done  by  various  investigators  on  the  causes 
of  these  unsatisfactory  lateral  oscillations.  Aerodynamic  balance  of  the 
rudder  is  a  probable  major  cause.  Reference  2  traced  the  poor  damping  qualities 
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of  its  test  airplane  (Navy  AD)  to  the  rudder  hinge  moment  characteristics. 
Reference  3  points  out  as  a  possible  cause  the  effect  of  fuel  sloshing  on 
airplane  dynamics.  Non-linearity  of  the  directional  stability  and  damping- 
in-yaw  derivatives  is  shown  to  be  a  possible  cause  in  reference  4.  Random 
separation  phenomena  such  as  flow  at  the  tail-fuselage  juncture  and  at  the 
tail  pipe  have  been  demonstrated  to  have  important  effects.  Reference  5 
shows  that  positive  inclination  of  the  principal  axis  can  have  considerable 
stabilizing  influence.  Here  wing  incidence  has  an  important  effect.   It  has 
been  demonstrated  theoretically  that  atmospheric  turbulence  can  initiate  and 
maintain  lateral  oscillations. 

The  designer  is  faced  with  an  imposing  problem  in  selecting  the  proper 
hinge  moment  parameters,  which  will  satisfy  control  force  and  other  require- 
ments, and  still  not  invite  snaking.  Reference  2  discusses  this  problem. 
Irreversible  power-boosted  control  systems  offer  a  promising  solution.  In 
any  case,  it  has  become  apparent  that  some  means  must  be  found  to  augment 
aerodynamic  dan^)ing.  The  increment  needed  is  considerably  in  excess  of  the 
amount  that  can  be  provided  by  an  increase  of  vertical  tail  area  or  tail 
length.  It  can  be  shown  that  the  tail  volume  coefficients  should  be  deter- 
mined by  the  requirements  for  static  sideslip  stability,  unsymmetrical  power, 
and  stalling  qualities,  and  that  other  methods  should  be  sought  to  add  to 
aerodynamic  damping. 

Although  hinge  moment  parameters  combined  with  friction  and  slack  in 
the  control  systems  explain  conventional  rudder  snaking  in  many  cases,  others 
are  on  record  in  which  the  rudder  has  been  locked  or  observed  not  to  move. 
Usually  the  calculated  and  tested  control  fixed  periods  and  times  to  damp 
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agree  quite  well.  The  case  under  consideration  in  this  paper  is  one  of  the 
latter  group.  The  manufacturer  felt  that  the  lateral  oscillations  were  trace- 
able to  non-linearities  in  the  derivatives  and  to  separation  effects.   It  was 
originally  thought  that  fuel  sloshing  was  a  contributory  cause,  but  baffling 
was  installed  and  found  to  have  no  noticeable  effect.   In  their  tests  they 
confirraed  the  effectiveness  of  a  "de-boost"  tab,  as  reported  in  reference  2« 
Since  there  are  no  known  aerodynamic  fixes  for  the  poor  damping  in  the  rudder 
fixed  condition,  an  automatic  stabilization  device  provides  the  most  promising 
solution. 

It  is  the  purpose  of  this  analysis  to  consider  the  effect  of  using  a 
servo  stabilizing  control,  and  to  select  an  optimum  setting  of  gain  and  spin 
axis  angle  for  this  application.  Two  systems  are  already  commercially  devel- 
oped (by  -Minneapolis-Honeywell  and  Boeing  for  the  XB-47,  reference  3,  and  by 
Lear  for  use  in  conjunction  with  the  F-5  autooilot,  reference  6).  The  theo- 
retical investigation  is  here  conducted  with  the  use  of  one  of  these  types  of 
servoraechanisros  in  mind. 

Reference  7  considers  stabilizing  devices  sensitive  to  various  quan- 
tities and  concludes  that  the  only  one  which  has  the  capability  of  satisfac- 
torily meeting  the  requirements  of  reference  1  is  an  autopilot  which  supplies 
rudder  control  proportional  yawing  velocity.  This  is  the  type  of  yaw  damper 
referred  to  above.  This  reference  eind  the  above  mentioned  systems  neglect 
the  problem  of  the  tilt  of  the  gyro  axis.  The  gyro,  which  is  the  rate  of  yaw 
sensing  device,  will  pick  up  components  of  roll  as  well  as  yaw,  and  at  only 
one  angle  of  attack  will  be  a  pure  yaw  sensing  instrument  as  is  customarily 
assumed.  This  fact  can  be  used  advantageously  to  effectively  gear  the  rudder 
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to  the  ailerons  by  making  the  yawing  moment  due  to  roll  positive.  This  elim- 
inates adverse  yaw  upon  entry  into  turns  and  rolling  maneuvers.  A  stabilizing 
device  can  be  designed  then  that  will  improve  the  damping  in  yaw  to  meet  the 
lateral  oscillation  requirements,  making  the  aircraft  spirally  stable  at  the 
same  time,  and  also  improve  the  response  to  the  rolling  control.  The  problem 
of  maintaining  equal  effectiveness  of  this  control  with  variations  in  speed 
is  also  discussed  in  this  analysis. 
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DESCRIPTION  OF  THE  AIRPLANE  AND  THE  DAMPER  SYSTEM 

The  airplane  under  consideration  is  a  single  seat,  midwing,  high 
performance  jet  propelled  fighter  with  a  conventional  control  system.  It 
has  a  two  piece  rudder  actuated  by  cables  fron  the  rudder  pedals,  not  power 
boosted.  The  physical  characteristics,  dimensions,  and  stability  derivatives 
were  supplied  by  the  manufacturer  and  checked  with  independent  wind  tunnel 
data  or  methods  of  calculation  as  specified  by  the  NAG A. 

The  damper  system  to  be  applied  in  this  case  would  be  known  as  a 
force  type  as  opposed  to  a  power- boosted  type.  This  means  that  the  servo 
drive  has  to  have  an  output  of  the  same  order  of  magnitude  as  a  human  pilot's 
control  effort.   In  the  other  case,  the  servo  drive  can  operate  into  an  arti- 
ficial feel  system,  boost,  or  irreversible  control  power  control  system.  Of 
course,  the  time  lag  of  this  force  type  system  is  not  increased  by  the  in- 
clusion of  the  time  lag  inherent  in  hydraulic  rudder  boost  systems.  Consid- 
erable range  in  power  size  is  available  to  meet  varying  requirements.  Time 
lags  for  a  force  type  system  can  be  reduced  to  less  than  one  tenth  of  a  second. 

The  proposed  damper   system  is  comprised  of  three  principal  parts, 
(l)  A  rate  gyro,  a  single  degree  of  freedom  gyro  mounted  with  its  spin  axis 
in  the  XZ  plane  of  the  aircraft  at  some  angle  to  the  body  reference  axis. 
It  is  constrained  by  a  spring  control  whose  restraint  is  caused  to  vary  to 
compensate  for  speed  changes.  The  gyro  produces  a  precession  displacement 
and  a  corresponding  electrical  signal  proportional  to  the  angular  rotation 
of  the  aircraft  in  space.  (2)  A  controller-amplifier  section,  which  receives 
the  rate  signals  from  the  rate  gyrot  amplifies  this  signal  and  imparts  phase 
lead.  It  contains  a  derivative  network  which  will  not  pass  steady  signals 
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from  the  rate  gyro.  This  network  prevents  the  damper  from  opposing  the 
pilot  in  a  steady  state  maneuver  such  as  a  turn  or  roll.  Any  transient  dis- 
turbances in  the  steady  state  maneuver  will  be  damped  out,  however,  and  the 
beneficial  damper  contribution  upon  entry  into  these  maneuvers  passes.  The 
gain  Eind  the  phase  load  are  adjustable  to  obtain  the  optimum  settings  for 
particular  applications.  The  compensating  spring  restraint  feature  could 
be  included  in  this  unit  instead  of  at  the  gyro.  A  typical  unit  would  also 
house  the  rate  gyro  in  its  chassis.  (3)  The  servo  drive  supplies  the  power 
to  move  the  rudder  upon  signal  from  the  eimplifier-contr oiler.   If  used  in  an 
irreversible  application  (which  it  is  not  in  this  case)  the  feedback  from 
the  load,  or  position  follow-up  would  not  be  required  since  the  artificial 
feel  spring  systems  will  recenter  the  load  when  the  signal  vanishes. 

The  total  weight  of  this  system  could  be  kept  to  within  fifteen 
pounds,  and  the  space  required  to  less  than  one  third  of  a  cubic  foot.  The 
gyro  should  be  located  near  the  center  of  gravity  so  as  to  escape  the  effects 
of  body  bending.   In  this  case  body  bending  should  be  a  negligible.  This 
problem  is  treated  in  reference  8. 

It  is  appropriate  that  the  distinction  between  an  autopilot  gind  a 
stabilizing  device  should  be  mentioned.  The  autopilot  has  authority  over 
the  aircraft,  presumably  by  the  choice  of  the  pilot.  He  may  overpower  it  or 
control  it  throu^  autopilot  control  settings,  but  the  autopilot  has  direct 
control  of  the  aircraft.  A  stabilizing  device,  such  as  the  proposed  yaw 
damper,  operates  in  conjunction  with  the  pilot*  the  pilot,  however,  has  con- 
trol and  is  superimposing  his  control  motions  over  the  deiraper.  The  damper 
can  be  thou^^t  of  in  the  same  sense  as  the  vertical  tail  in  that  it  is  a  part 
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of  the  aircraft  which  ia  in  continual  operation  contributing  stability.   It 
is  the  opinion  of  the  author  that  it  should  be  an  absolute  requirement  that 
a  stabilizer  shall  have  no  feedback  to  the  control  system  of  the  pilot  1  the 
pilot  should  not  be  conscious  of  the  damper's  operation  except  that  he  feels 
that  the  flying  qualities  are  improved,  and  that  he  is  able  to  fly  a  more 
hl^ly  damped  and  smoother  path.  This  isolation  can  be  easily  accomplished 
in  an  irreversible  control  system.  In   the  application  under  consideration 
it  poses  a  difficult  problem  that  can  probably  be  best  solved  by  dividing  the 
two  piece  rudder  between  the  pilot  and  the  damper.  That  is,  using  a  control 
surface  for  the  damper  alone.  No  solution  in  which  the  pilot  can  feel  the 
damper  in  operation  throu^  the  rudder  pedals  should  be  considered  satisfactory. 
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UTERAL  DYIIAMIC  CHARACTERISTICS 

The  first  step  in  the  problem  is  to  determine  what  the  calculated 
lateral-directional  oscillation  characteristics  of  the  aircraft  are. 

Three  oonfiguratione  of  the  aircraft  have  been  considered. 
I.  The  raaximuin  cruise  configuration,  at  35,000  feet,  clean,  at  a  l^ch 

number  of  .68. 
II.  The  505&  Normal  Rated  Power  configuration  at  10,000  feet,  clean,  at 

a  Mach  number  of  .32. 
III.  The  Povrer  Approach  condition  with  flaps  and  gear  down  and  with  power 

for  level  flight  at  a  reasonable  approach  speed. 

These  three  conditions  give  a  sufficient  range  of  altitude  and  speed 
to  reasonably  judge  the  lateral  dynamics  of  the  aircraft  with  and  without  a 
damper  system.  Configuration  II  was  selected  because  pilots  reports  have 
indicated  that  at  that  particular  altitude,  at  ^iach  numbers  below  .6,  the 
Dutch  Roll  or  snaking  has  been  particularly  noticeable.  Conditions  I  and 
III  represent  the  extreme  opposite  ends  of  the  speed  and  altitude  range  in 
configurations  of  operational  interest. 

The  equations  of  motion  as  presented  in  reference  9  have  been  used, 
with  the  exception  that  the  effects  of  the  product  of  inertia  have  been  in- 
cluded only  in  condition  I.   In  this  condition  "^X     is  negative..  In  the  other 
two  conditions  T[      is  positive  and  it  is  conservative  to  neglect  the  product 
of  inertia.  For  this  particular  aircraft  the  product  of  inertia  was  found 
to  have  little  effect.  These  calculations  are  reviewed  in  the  section  on 
calculations.  The  data  is  presented  in  Table  II  and  plotted  in  figure  2. 

All  the  points  lie  in  the  unsatisfactory  region.  Flight  tests  have 
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shovm  good  correlation  with  these  points.  The  calculations  are  control 
fixed,  indicating  an  unsatisfactory  Dutch  Roll  characteristic.  The  boundary 
between  the  satisfactory  and  unsatisfactory  regions  of  figure  1  of  reference 
1  is  currently  under  investigation.  The  calculated  points  show  the  aircraft 
to  be  in  a  region  which  has  been  called  tolerable,  meaning  tiring  or  un- 
pleasant, but  not  necessarily  dangerous.  The  calculations  show  that  a  con- 
dition exists  which  would  ceirtainly  be  einnoying  and  detrimental  to  the 
carrying  out  of  the  airplane's  mission. 

Whether  this  condition  is  of  sufficient  severity  to  warrant  a  change 
or  modification  in  already  operational  units  is  doubtful.  That  is  a  question 
of  opinion  of  the  pilots  operating  the  airplane.  /U.1  indications  are  that 
future  designs  are  going  to  need  automatic  stabilization!  analyses  of  this 
type  will  be  called  for.  Applications  to  longitudinal  motions  are  also  in 
use.  Future  airplanes  may  well  use  automatic  stabilization  both  laterally 
and  longitudinally  with  the  pilot  imposing  his  control  efforts  over  the 
automatic  controls. 
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CONTRIBUTION  OF  TEE  YiW  DAKiPER 

Referririi^  to  figure  lb  the  relation  between  the  body  reference 
axis,  the  gyro  spin  axis,   and  the  wind  or  stability  axis  is  shown.     Sub- 
scripts R  (reference),   G  (gyro),   and  V  (stability)  are  used  in  the  following 
development. 

For  simplification  the  time  lag  of  the  yaw  damper  system  is  first 
assumed  to  be  zero.     It  will  be  later  shown  that  the  critical  time  lag  of 
the  aircraft-stabilizer  system  is  considerably  hi^er  than  the  time  lag  of 
any  mechanism  under  consideration. 

The  equations  of  motion  are  written  with  respect  to  the  wind  or 
stability  axes.     The  derivatives  contributed  by  the  damper  must  therefore 
be  related  to  the  V  axis  for  inclusion  in  the  equations.     The  rudder,  how- 
ever,  introduces  yawing  and  rolling  moments  about  the  reference  or  body  axis 
R.     Also  the  yawing  and  rolling  moments  so  introduced  are  proportional  to 
the  rate  of  yaw  as  measured  with  respect  to  the  gyro  axis  G. 

Then  the  development  of  the  contributions  of  the  damper  is  as  follows: 

^  HdT^!       ^^^^rfid/sec     °°"*^°^  gearing  or  gain. 
G 

^  n  B  ^£       ^  %  ^ 

r  r 

r  r 


but 


and 


H^Q-        4^^003?        +^Slllf 
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Cj2  *  (Ci?)y  =  (C^  )q  cos  a  +  (C^)q  sin  a 
substituting  into  the  above  equations f 

C     «KC         [H^^cosf     +^  sin  |     ]  cos  a 

+K  C^       C  4^  y  cos  f     +  0^  sin  I     ]  sin  a 

Og^  ■»  Yl  C^       [  ^  y  COS  I     +  0^  sin  I     ]  COS  a 
r 

+K  C^      [  4^  Y  cos  I     +  ^  sin  I     ]  sin  a 

These  expressions  can  be  seen  to  have  the  form  of  total  derivatives,  where 

dC  .  dC         . 

^^n'    df    ^^     ^     df  ^ 

dC  o         .  dCjj        , 


therefore 


AC         «  K  cos  t      [C       +     C  <,       a] 


AC        -Ksinf     [C      +    Cj^    a] 
'^  r  r 


AC/,.     »  K  coa  I     [C.     +  C         a] 

r  r 

AC,.     '  K  sin  I      [C(7     +  C         a] 
'^                                  r  r 


since 
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C         =     C         :^     '    C       r*       etc. 
2V 


2V  b 

AC     =    —■  K  cos  I     C 


AC     =    §^  K  3in  ^     C 
n_        b  ^       n< 

r 


P 


neglecting  the  product     o.  Ca       ,  which  is  small  compared  to     C 

0  6 

r  r 

The  daraper  contributions  are  seen  to  be  directly  proportional  to 
true  airspeed.  The  rudder  power  also  is  seen  to  vary  with  CLLtitude  from 
Table  II. 

The  rudder  power  coefficient  has  a  negative  sign,  and  the  rolling 

coefficient  due  to  rudder  deflection  has  a  positive  sign.  It  can  be  seen, 

therefore  that  if  f   has  a  negative  sign  the  increnaents  AC   and  AC^,  will 

'  ^  ^^  ^ 

be  positive,  and  the  increments  AC    and  AC^   negative.  These  are  the 
desirable  signs  for  these  contributions  by  the  damper.   It  is  therefore 
concluded  that  the  gyro  axis  should  be  so  set  that  it  is  always  below  the 
stability  axis,  (  f   a  negative  angle)  eis  shown  in  figure  Ic. 

The  increments  /!^n       and  ££/,  .   have  little  effect  on  the  lateral 
dynamic  calculations,  and  sum  [Cn   +  C   a]  is  considerably  less  than 

T  T 

[C   +  Cn  a] I  for  these  reasons  the  conservative  assumption  is  made  in 

r    r 

further  calculations  that  only  AC     and  AC  ,  are  supplied  by  the  yaw 

damper. 

The  increments  added  by  the  yaw  damper  are  then  assumed  to  be 
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AC         »    K  cos  I     C 
AC.     '    K  sin  I      C 

The  optimum  settings  of  K  and  |  =  f(a,D)  are  discussed  in  a 

later  section. 

2V 
Reference  10  treats  the  effect  of  C  {-r   C  )  on  lateral  dynamic 

^p  ^  '  ^^ 
stability,  llaking  C  ,  positive  was  shown  to  have  a  beneficial  effect  in 

^^ 

shifting  the  points  on  a  plot  such  as  figure  4  to  the  satisfactory  area. 

As  K  and  -f   are  increased,  the  positive  increment  of  C  .  lengthens  the 
period  at  substantially  the  same  time  to  damp  to  one-half  amplitude,  improv- 
ing the  cycles  to  damp  to  one-half  amplitude.  This  increment  in  +C   must 
be  limited,  however,  since  too  great  an  amount  leads  to  a  long  period  un- 
stable oscillation  of  the  roots  that  were  identified  as  the  rolling  converg- 
ence and  the  spiral  mode.  These  roots  first  become  a  long  period  stable 
oscillation  and  then  as  C   is  further  increased  the  oscillation  becomes 
unstable.  This  shows  that  the  product  VK  sin  f  must  be  limited  so  as  to  be 
satisfactory  in  configurations  as  contrasting  as  I  and  III, 

The  reasons  for  seeking  +C  ,  contributions  from  the  damper  are  to 
assist  in  correcting  the  Dutch  Roll  oscillations  and  to  oppose  the  adverse 
yaw  introduced  by  the  ailerons  and  the  normal  negative  derivative  (adverse) 
of  yaw  due  to  roll.  By  overcoming  the  adverse  yaw  in  this  manner,  entry 
into  rolls  and  turning  fli^t  can  be  considerably  improved.  This  benefit 
can  only  be  enjoyed  on  the  entry  into  these  maneuvers,  since  the  steady  yawing 
velocity  signal  will  be  intentionally  blocked  to  keep  the  yaw  damper  from 
opposing  a  steady  turn  initiated  by  the  pilot. 
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C   is  seen  to  have  its  greatest  influence  in  the  b,  term  of  the 
stability  quart ic.  As  this  derivative  becomes  more  positive  it  drives  b, 
to  zero.  As  this  occurs  the  rolling  and  spiral  modes  first  become  a  stable 
oscillation,  and  finally  a  long  period  unstable  oscillation.  The  limiting 
value  of  - 1   for  the  fixed  K  was  then  calculated  by  setting  b,  to 

zero,  keeping  all  derivatives  except  C   and  C  ,  the  same,  and  introduc- 

^    n^      n^. 

Ing  these  as  f ( |  ).  The  limiting  P  determined  in  condition  I  is  shown 

in  figure  3.  Since  a  varies  over  a  small  known  range,  and  5  «  a  -  P  , 

it  is  possible  to  estimate  the  range  of  effective  f   .  See  figures  Ic,  3. 

££,       emd  tC^    ,  the  non-dimensional  forms  of  the  damper  contribu- 
n       n 
r        p 

tions,  were  seen  to  vary  with  speed.  Therefore  in  order  to  maintain  equal 
effectiveness  of  the  stabilizing  control  throu^out  the  speed  range,  K 

must  be  made  to  vary  inversely  with  speed.  One  convenient  way  to  do  this 

1/2 
is  to  apply  the  inverse  of  q  '   to  the  spring  restraint  of  the  gyro, 

thereby  modifying  the  electrical  signal  from  it.  The  V's  will  cancel,  but 

-1/2 
a  varying  factor  of  cr  '   remains  that  will  have  to  be  calibrated  out  as 

well  as  possible.  A  close  approximation  to  the  desirable  variations  ip  K 
with  the  three  conditions  is  shown  in  Table  III  (K  to  make  C.  /g*  .5).  Com- 
paring these  values  with  those  shown  in  Table  II  (calculated  by  using  a 

-1/2 
q  '   variation)  it  is  seen  that  the  contributions  tend  to  be  too  large  at 

altitude.  Figures  2,  3,  and  4  confirm  this.  Applying  dynamic  pressure  from 

a  pitot  source  seems  to  be  the  most  convenient  way  of  compensating,  however. 

A  non-linear  spring  and  bellows  with  a  pitot  tube  source  such  as  are  used 

in  artificial  feel  systems  would  probably  be  the  best  solution.  This  penalty 

in  gadgetry  must  be  paid  to  realize  the  benefits  of  the  positive  yaw  due  to 

roll  introduced  by  the  damper.  Positive  C   must  be  limited,  and  the  goal 

P 
should  be  similar  performance  at  all  altitudes  and  speeds. 
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SETTING  GIRO  AXIS  ANGI£  AND  GAIN 

The  question  of  just  how  damped,  a  motion  the  pilot  would  like  to 
have  is  still  under  considerable  discussion.     There  is  no  question  that  in 
a  system  of  this  sort  it  would  be  possible  to  add  enou^  damping  with  a 
stabilizer  to  the  aerodynamic  damping  to  critically  damp  all  lateral  oscil- 
lations.    However,  the  standards  for  the  performance  and  stability  of  the 
aircraft-stabilizer  servomechanism  are  not  as  yet  established.     See  refer- 
ences 11  and  12, 

This  analysis  was  made  on  the  assumption  that  the  gain  should  be 
the  minimum  required  to  bring  the  cycles  to  one-half  amplitude  to  one  half. 
There  is  evidence  to  show  that  the  motion  will  then  appear  critically  damped 
to  the  pilot.     Reference  8. 

A  Routh's  discriminate  investigation  using    K     and    |      as  variables 
is  possible,  but  exceedingly  tedious.     Determination  of  the  optimum  values 
of  these  variables  by  a  plot  of     C,  /„     vs    0     with     K     as  a  parameter  is 
shown  in  figure  3.     This  approach  is  relatively  easy  euid  equally  informative. 

Referring  to  figure  3,   curves  for     K  «  1     and     K  «  2     in  condition 
III  are  shown.     The  higher  value  of     K     does  not  give  sufficiently  better 
results  to  warrant  the  hi^er  gain  settings.      It  is  in  this  condition  (sea 
level  at  approximately  1.20  V   .    .^)  that  the  basic  value  of  the  gain,   or 
control  gearing,   is  established.     At  hi^er  speeds   (or  q's)  the  gain  is 
reduced  as  already  explained. 

Selecting  the  basic  value  of     K  =  1     in  condition  III,  the  effective 
value  in  condition  I  becomes   .5,   in  condition  II,    .817.     The  curve  for  con- 
dition I  for  this  effective  K  is  shown  up  to  the  limiting  value  of    [P     , 
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determined  as  discussed  in  the  previous  section.  In  this  particular  case 
there  is  no  limit  to  the  gyro  angle  in  condition  III.  In  other  words  with 
gains  below  2,  too  much  +C  cannot  be  introduced.  The  effect  of  holding 
K  constant  and  inclining  the  gyro  axis  to  the  point  that  damper  contribu- 
tion is  entirely  ££,  is  shown  in  figure  4.  The  points  of  |  =  0  and 
f  "  -90  on  each  curve  contrast  the  effects  of  stabilizers  sensitive  to 
pure  yaw  or  pure  roll. 

Selecting  a  gyro  angle  of  30  and  K  «  1,  the  range  of  values  of 
C-  /„  for  a  variation  of  angle  of  attack  of  fifteen  degrees  is  shown.  This 
setting  of  D  and  K  is  taken  for  further  investigation  since  it  is  the 
maximum  practical  gyro  angle  at  this  K  that  will  not  introduce  an  unstable 

oscillation  due  to  too  much  +C   in  condition  I.  Some  lower  angle  might 

n  w    u 

P 

prove  more  desirable  from  fll^t  tests.  This  is  a  question  of  how  much 

positive  yaw  due  to  roll  is  best  for  entry  into  rolls  or  turns.  From  figures 

8,  9,  and  10  it  can  be  seen  that  the  response  to  the  aileron  control  is 

raidically  modified. 

In  general,  a  quick  procedure  for  selecting  gyro  angle  and  gain  is 

as  follows  I 

(1)  Using  the  single  degree  of  freedom  equation  in  yaw,  a  close 
estimate  of  the  total  amount  of  damping  required  to  damp  to  any  given  C.  /_ 
is  easily  obtainable.  From  this,  the  amount  of  AC  ,  that  is  to  be  added 
by  the  damper  can  be  estimated.  Knowing  the  rudder  power,  the  approximate 
gain  to  get  this  increment  of  damping  can  be  found, 

(2)  From  the  adverse  yaw  and  yaw  due  to  roll  coefficients  the  total 
amount  of  ^t£.   ^   desired  can  be  estimated.  Using  the  gain  calculated  above. 
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the  required  I   to  get  the  increment  in  +C   estimated  above  can  be  solved 
for,  using  the  expression  obtained  under  the  daaiper  contribution  section. 
For  the  gyro  angle,   l^  «  a  -  f  . 

The  compensating  device  on  the  gyro  will  keep  these  increments  essen- 
tially constant  with  speed  and  altitude.  To  check  to  see  if  an  instability 
has  been  introduced,  select  a  high  speed  condition  where  the  final  ratio  of 
C  .  to  C  .   is  high  (usually  at  high  altitude  and  speed)  and  calculate  the 
b,  term  of  the  stability  quart ic. 

Final  checks  using  the  three  degree  of  freedom  equations  and  involv- 
ing a  separate  solution  of  the  quartic  for  each  point  as  was  done  in  figure 
3  can  then  be  made  for  a  critical  configuration.  The  assumption  of  zero 
time  lag  will  have  to  be  kept  in  mind  and  later  checked.  At  this  point  it 
probably  would  be  desirable  to  begin  fli^^t  testing  the  system  for  the  final 
selection  of  these  paramenters. 

In  order  to  prevent  any  feedback  through  the  pilots  rudder  control 
system,  the  yaw  damper  for  this  application  will  probably  have  to  have  its 
own  rudder  area.  The  calculations  in  this  report  have  been  made  on  the 
assumption  that  the  whole  rudder  was  available  for  the  yaw  damper.  This 
can  be  done  when  the  movements  of  the  rudder  caused  by  the  danger  can  be 
masked  from  the  pilot  through  a  boosted  or  artificial  feel  control  system. 

c   =  -c   ^ 

»5,  \,  ^^ 


and 


Sph     ^ 
r       r  r 


j 
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The  terras   involving    C/j     have  been  neglected  in  the  develoninent 

r 
of  the  contribution  of  the  yaw  damper.     Since 

AC         =     K  cos   I      C 

r 

AC         «     K  ain   I      C 

^  r 

the  same  results  can  be  obtained  by  keeping  the  product  of  KC     constant: 

r 

or  as  the  area  Sj.,  for  the  dsunper,  is  reduced,  the  gain  K  is  increased. 

Note  that  the  vertical  location  of  the  area  for  the  damper  is  a  negligible 
effect.  The  maximum  practical  value  of  K  obtainable  with  this  type  of  in- 
stallation (probably  2  or  3)  will  limit  the  mininum  size  of  the  damper  rudder 
area. 

In  surnmary,  the  optimum  gyro  angle  is  chosen  as  30  degrees,  and  the 
optiraun  gain  as  1  at  sea  level  and  approach  speed.  This  setting  has  the 
following  results: 

(1)  Moves  the  points  shown  in  figure  2  from  the  unsatisfactory 
region  to  the  satisfactory  damping  region  where  C.  /„  is  less  than  .5. 

(2)  Adds  a  sufficient  amount  of  AC  .  in  each  case  to  make  the 
aircraft  spirally  stable  in  all  conditions. 

(3)  Adds  a  +C   increment  of  an  approximately  constant  value  to 
all  conditions,  which  is  considerably  larger  than   the  total  adverse  yaw  due 
to  roll  and  the  rolling  control. 

The  beneficial  effect  of  the  additional  damping  is  shown  in  the 
response  data  plotted  in  figures  5,  6,  and  7.  The  improved  response  to  the 
rolling  control  is  shown  in  the  response  data  plotted  in  figures  3,  9,  and 
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10. 

The  effect  of  varying  P     at  a  constant  K  for  condition  III  is 
shown  in  figure  4. 

Figure  3  shows  the  selected  points  in  the  three  conditions,  and  a 
C,  /p  vs  U>  plot  for  constant  K  for  conditions  I  and  III  including  the 
range  of  angles  of  attack. 

All  calculations  involving  the  contributions  from  the  yaw  damper 
were  made  without  using  the  small  angle  assumptions  or  neglecting  the  influ- 
ence of  roll  due  to  rudder  deflection.  The  simplifying  assumptions  made  in 
the  analysis  of  the  damper  contributions  were  found  to  be  justified. 


I 
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ADV/iNTACE  OF  YAW  DAMPER  STABILIZATION 

It  was  pointed  out  in  the  Introduction  that  aerodynamic  damping 
should  be  augmented  by  some  other  means  than  increasing  the  tail  length  or 
area.  Increased  tail  length  is  limited  by  operational  factors  such  as  car- 
rier elevator  sizes  for  naval  aircraft.  The  simount  of  vertical  tail  area 
also  has  a  practical  lirdtation  which  precludes  the  increases  in  damping 
that  are  desired. 

The  yaw  damper  has  the  advantage  of  increasing  C   without  affect- 

r 
in^  the  value  of  C  .  Tail  size  can  then  be  chosen  without  consideration 

of  the  requirement  of  hi^  damping,  and  C   can  be  set  from  consideration 
of  static  stability,  stalling  qualities,  etc.  Dismissing  the  possibility 
of  adding  to  conventional  tail  length,  the  vertical  tail  area  would  have  to 
be  increased  by  a  factor  of  about  four  to  achieve  the  same  results  in  damp- 
ing that  the  addition  of  this  stabilizer  will. 

Reference  8  points  out  that  with  the  XB-47  it  was  physically  impos- 
sible to  change  the  effective  dihedral  negatively  enou^  to  correct  the 
lateral  dynamics  because  of  ground  clearance,  and  that  the  tail  area  would 
have  to  be  at  least  doubled  to  add  the  necessary  damping.   It  was  explained 
that  in  general,  while  increasing  tail  size  increases  damping,  the  gain  in 
smooth  flying  qualities  is  questionable.  While  the  tail  size  increase  adds 
to  the  damping,  the  excitation  of  the  aircraft  due  to  side  gusts  increases 
in  the  same  proportion. 

This  XB-47  report  uses  the  servomechanisms  approach  to  stability 
analysis  as  outlined  in  references  11  and  12.   In  one  polar  stability  plot 
White  shows  that  the  effect  of  increasing  C   and  C   together  (as  in  the 


I 
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case  of  increasing  tail  area)  results   in  little  or  no  improvement  in  damping 
because  the  resonant  frequency  of  the  aircraft  has  been  increased. 
Consider  the  single  degree  of  freedom  equation  in  yaw, 

*^  r 

This  equation  may  be  transformed  into  the  standard  servomechanisms 
form  using  the  following  definitions: 

-C  . 

damping  ratio  «3  «  ^j  ^ 

z  o 

C 

undamped  natural  frequency    co  «   ,  ^  ' 

z 


rudder  effectiveness  K 


C 


r    c 
n 


P 


and  the 


The  equation  of  motion  then  becomes 

•^  +  2%  CO  4^  +  w^  4^  '    K,  0?  S: 

o              o  •  ■  o 

damped  natural  frequency  u)  *  ^n/^  "  3 


^yi  -  25 


resonant  frequency  ui  *=  co 

For  the  usual  values  of  5   ,  o)  ,  co  ,  and  O)   may  be  considered 

*  o   r      n  ^ 

equal.  White,  reference  8,  considering  the  single  degree  of  freedom  case, 
sets  up  his  criterion  of  damping  as  the  ratio  of  dynamic  yaw  response  occur- 
ring at  u)  ,  to  the  effective  static  yaw  response  considering  a>  *  a)  . 
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The  measure  of  the  damping  then  becomes   inversely  proportional  to  the  square 

of  the  xmdamped  natural  frequency,   hence   inversely  proportional  to     C        , 

2  2/  P 

since    oj    «i    o)      =    C       /I     . 
r  o  n^  '    z 

P 

Also  from  the  single  degree  of  freedom  equation, 

r    ,  «695  ^1  «  '6  ^ 
^1/2     2tr    5 

White  states  that  within  the  limits  of  engineering  accuracy,  the 
single  degree  of  freedom  equations  may  be  used  for  calculations  involving 
the  design  of  the  yaw  damper  and  to  obtain  C^  /„  • 

Results  from  single  degree  of  freedom  calculations  are  included  in 
Table  III. 


PRINCETON  UNIVERSITY  PAGE   23 

AERONAUTICAL  ENGINEERING  LABORATORY  report  172 


TRANSIENT  RESPONSE  INVESTIGATIONS 

In   order  to  graphically  present  the  contrast  between  the  transient 
response  of  the  aircraft  with  and  without  the  damper  system  operating,  the 
following  calculations  were  made.  For  condition  II  the  responses    vP  « 
P/P  ,  emd  0/p   were  found  for  a  disturbance  in  sideslip  effecting  each 
of  the  three  equations  of  motion.  The  function  introduced  is  a  unit  step 
function  in  sideslip  applied  as  a  step  function  at  t  =  0,  and  neutralized 
as  a  step  function  at  t  "  1.  Condition  II  was  selected  because  in  this 
configuration  pilots  have  frequently  reported  unpleasant  lateral  oscilla- 
tions. The  forcing  function  used  is  an  attempt  to  simulate  a  disturbance 
that  might  be  initiated  by  turbulence  or  by  a  quick  displacement  in  sideslip 
however  introduced. 

The  airplane  in  the  original  condition  is  seen  to  have  persistent 
oscillations  in  all  three  plots  (figures  5,  6,  and  7).  The  period  is  long 
enough  so  that  the  average  pilot  could  damp  them  out  with  the  proper  appli- 
cation of  controls;  however,  requiring  the  pilot  to  continually  oppose  these 
motions  is  tiring  and  detrimental  to  carrying  out  his  mission.  There  is  a 
considerable  improvement  in  yawing  velocity  and  sideslip,  but  the  rolling 
characteristic  with  the  damper  in  operation  should  be  investigated  in  flight 
test.  Too  high  a  ratio  of  roll  to  sideslip  is  particularly  unpleasant  to 
the  pilot.  From  data  available  in  other  applications,  oscillations  that  damp 
to  half  amplitude  in  less  than  half  a  cycle  appear  dead  beat  to  the  pilot. 

From  this  data  and  the  rudder  hinge  moment  parameters,  the  rudder 
pedal  force  that  must  be  supplied  by  the  servo  power  section  of  the  damper 
can  be  estimated.  An  estimate  of  the  amount  of  rudder  throw  required  of  the 
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damper  rudder  surface  can  also  be  made.  The  final  setting  of  the  stops 
limiting  the  throw  of  the  damper  rudder  surface  must  be  determined  experi- 
mentsdly  since  the  control  required  will  depend  upon  the  characteristics  of 
the  servo  system  as  well  as  the  airplane  responses  and  the  firequency  range 
to  be  deiraped. 

The  response  to  an  aileron  step  function  is  presented  in  figures 
8,  9,  and  10.  Condition  III  is  investigated  since  it  is  probably  the  most 
unfavorable.  Adverse  yaw  has  been  included  in  the  calculations.  The  large 
positive  yaw  due  to  roll  introduced  by  the  yaw  damper  is  seen  to  drastically 
change  the  initial  response  upon  entry  into  the  rolling  maneuver.  The  air- 
plane no  longer  starts  to  yaw  in  the  wrong  direction  for  the  turn,  and  the 
roll  builds  up  rapidly.  The  sideslip  builds  up  at  a  slightly  faster  rate 
also,  but  not  rapidly  enou/^  to  offset  the  good  effects  in  yaw  and  roll. 
The  derivative  network  in  the  yaw  damper  circuit  will  have  an  RC  circuit 
with  a  time  constant  so  designed  as  to  only  pass  the  first  portion  of  a  second 
of  the  contribution  of  the  damper  so  that  it  will  have  the  effect  of  assist- 
ing in  the  rapid  build-up  of  roll  without  the  adverse  initial  yaw,  and  still 
not  oppose  the  steady-state  maneuver.   In  this  manner  the  damper  should  con- 
siderably improve  the  response  to  the  rolling  control. 

From  figure  5,  (the  response  in  yawing  velocity  to  the  disturbance 
in  sideslip)  the  maximum  yawing  velocity  per  degree  disturbance  is  seen  to 
be  approximately  two  degrees  per  second  per  degree  sideslip.  Then  ten  degrees 
of  rudder  would  be  required  for  a  gain  of  one  and  an  initial  disturbance  of 


five  degrees,  since  S      *  jd^  /d  4^  H^ . 

Considering^  a  ten  degree  throw  in  either  direction  as  the  limit  for 


I 
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damper  applied  rudder  deflection,  it  can  be  seen  that  the  "\/ith  damper" 
curves  shovm  in  the  aileron  response  plots  will  be  cut  off  by  the  rudder 
reaching  its  stops.  This  cut-off  point  will  be  determined  by  the  amount  of 
aileron  applied  since  figure  9  is  presented  per  degree  deflection.  When 
the  product  of  the  gain  and  the  yawing  velocity  exceeds  the  available  full 
deflection,  the  rudder  will  be  initially  held  against  the  stop  with  a  ten 
degree  favorable  deflection.  The  RC  circuit  will  then  cut  out  the  damper 
rudder  contribution  as  the  steady  signal  is  sensed  by  the  gyro.  Should 
there  be  a  tendency  to  oscillate,  the  aervomechanism  will  damp  the  variations 
about  the  steady  maneuver  out. 

The  curves  showing  the  improved  performance  with  the  damper  oper- 
ating (figures  3,  9,  and  10)  will  then  be  modified  by  servo  rudder  stop 
setting.  Startling  improvement  in  the  entry  into  rolling  maneuvers  should 
still  be  obtainable  with  limited  rudder  throw. 

Two  factors  should  be  here  reiterated.  First,  these  calculations 
were  made  for  the  maximum  value  of  gyro  axis  inclination,  and  hence  positive 

C   increment,  possible  at  the  selected  value  of  gain.  In  other  words  these 

P 
aileron  responses  are  for  the  maximum  gearing  of  the  rudder  to  the  ailerons. 

Some  smaller  angle  of  gyro  axis  inclination  mi^t  prove  more  desirable  from 

fli^t  test  or  more  exhaustive  calculations.  Second,  the  device  introduced 

f 

to  assure  nearly  constant  contributions  from  the  damper  with  variations  in 
speed  and  altitude  will  automatically  reduce  the  gain  at  hi^er  velocities. 
This  should  prevent  destructive  rudder  motions  as  it  hits  its  stop  after  the 
gyro  senses  a  large  rolling  velocity. 


I 
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DETERMINATION  OF  CRITICAL  Tlf^  UQ 

In  all  the  calculations  to  this  point  the  assumption  has  been 
made  that  the  time  lag  of  the  damper  system  is  zero.  The  critical  time  lag 
of  the  aircraft-autopilot  combination  is  now  investigated  to  determine  how 
much  of  a  limitation  this  assumption  imposes. 

The  frequency  response  method  outlined  in  reference  14  is  applied. 
The  critical  time  lag  of  the  system  is  defined  aa  the  time  lag  that  results 
in  a  neutrally  stable  or  steady-state  oscillation.  The  steady-state  response 
K.  «  j'^/^  I  o^  ^^®  airplane  to  a  sinusoidal  forcing  function  of  unit  ampli- 
tude was  first  calculated.  The  calculations  were  made  for  condition  III 
since  in  this  condition  the  ratio  of  the  damping  to  the  critical  damping 
was  lowest.  The  gain  of  the  yaw  damper  is  knownt  having  been  set  at  K«l 
for  this  condition.  The  gain,  K.  is  the  amplitude  of  the  control-surface 
oscillation  produced  by  the  yaw  damper  mechanism  in  response  to  the  oscil- 
lation of  the  airplane  yawing  velocity.  Intersecting  K  =  j l/K  j   with  the 
airplane  gain  or  amplitude  plot  the  maximum  permissible  time  lag  in  the 
system  can  be  determined.  The  slope  of  the  line  throu^  the  point  where  the 
phase  angles  are  equal  is  a  measure  of  this  critical  time  lag.  The  condi- 
tion for  stability  is  that  at  each  angular  frequency  where  the  phase  angles 
are  equals  the  ratio  of  the  amplitudes  of  the  airplane  to  the  autopilot  must 
be  less  than  one. 

From  figure  11  the  critical  time  lag  of  the  system  chosen  is  seen 
to  be  1.2  seconds.  Critical  time  lags  for  gains  of  two  and  three  are  shown 
also.  For  a  damper  gain  of  three,  the  limiting  time  lag  is  .7  considerably 
in  excess  of  that  inherent  in  typical  force  type  systems.  The  critical  lag 


I 
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shovm  in  reference  3  is  a  little  over  .5  seconds  for     K  «  1.     For  servo- 
mechanisms  with  larger  time  lags   introduced  by  boost  systems,   this  calcu- 
lation would  undoubtedly  show  a  closer  margin.      In  this  case,  the  time  lag 
of  the  system  is  not  a  limiting  consideration. 


' 


t 
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CONCLUSIONS 

A.a  a  result  of  the  theoretical  investigation  of  the  application 
of  a  yaw  damping  aervomechanism  to  the  aircraft  under  consideration,  the 
following  conclusions  are  reached; 

(1)  This  plane  does  not  satisfactorily  meet  the  requirements  of 
reference  1  in  regard  to  lateral-directional  oscillations.  Control-fixed 
calculations  predict  that  the  airplane  will  exhibit  a  "Dutch  Roll"  in  all 
the  conditions  investigated  to  a  degree  that  could  be  classified  as  objec- 
tionable but  not  dangerous.  These  calculations  substantiate  pilots  reports 
of  "snaking"  or  "Dutch  Roll." 

(2)  By  installing  the  yaw  damper  system  as  proposed,  the  aircraft 
will  then  meet  the  requirements  of  reference  1.  The  lateral  oscillations 
will  appear  nearly  critically  damped  to  the  pilot,  and  he  will  be  able  to 
carry  out  the  mission  of  the  aircraft  more  effectively  with  less  effort. 

(3)  The  damper  will  introduce  enou^  damping  in  yaw  to  make  the 
aircraft  always  spirally  stable.   It  will  also  introduce  a  constant  positive 
yawing  moment  due  to  roll  which  will  improve  the  turning  and  rolling  per- 
formance. The  system  will  be  so  designed  as  to  improve  entry  into  steady 
turning  or  rolling  maneuvers,  but  not  oppose  them. 

(4)  Damping  introduced  by  a  servoraechanism  which  applies  rudder 
control  proportional  to  yawing  velocity  is  probably  the  most  promising  if 
not  the  only  satisfactory  way  to  augment  the  aerodynamic  deunping.  This  sta- 
bilizer introduces  damping  without  increasing  the  resonant  frequency  (or 
the  effective  static  yaw  stability),  which  insures  results  not  obtainable 

by  increasing  vertical  tail  area. 


I 
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(5)  For  this  particular  installation  the  optimum  settings  of  gyro 
spin  axis  angle  and  damper  gain  are  U>  «  30  and  K  =  1.  These  values  have 
been  checked  in  a  limiting  condition  and  the  amount  of  positive  yaw  due  to 
roll  introduced  is  within  the  limit  that  would  cause  a  long  period  unstable 
oscillation.  Some  lower  value  of  U  that  would  still  produce  a  +C  ,  might 
prove  more  desirable  from  flight  tests  of  the  effect  on  rolling  maneuvers. 
In  all  cases,  the  spin  axis  of  the  gyro  should  be  kept  below  the  stability 
or  wind  axis.  Figure  Ic. 

(8)  It  should  be  required  in  the  case  of  a  stabilizing  system 
(as  opposed  to  an  autopilot)  that  no  feedback  to  the  pilots  controls  be 
allowed.   In  this  case  the  most  practical  solution  appears  to  be  the  use  of 
a  separate  rudder  area  for  the  yaw  damper  mechanism,  probably  splitting  the 
two  piece  rudder  between  the  pilot  and  the  damper.   A  differential  linkage 
or  gearing,  with  the  pilot  and  stabilizer  aia  independent  inputs  with  pro- 
vision for  balancing  out  the  servo  applied  hinge  moments  could  be  developed. 

(7)  Calculations  show  that  the  damper  system  as  proposed  will 
improve  the  transient  response  to  a  disturbance  in  sideslip,  and  also  improve 
the  response  to  an  aileron  step  function  input. 

(3)  The  critical  time  lag  of  the  stabilizer-airplane  system  is 
considerably  in  excess  of  the  maximum  time  lag  to  be  expected  of  any  servo- 
mechanism  of  the  type  oroposed.  The  assumption  of  zero  time  lag  in  the  cal- 
culations proved  to  be  reasonable. 

White,  in  reference  3, has  an  excellent  presentation  for  the  invest- 
igation of  the  stability  of  the  stabilizer-airplane  combination  which  he  calls 
the  inverse  frequency  response  method.  It  is  felt  that  this  type  of  approach 
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should  be  used  after  the  experimental  frequency  response  data  of  the  actual 
servomechanisn  to  be  used  is  available.  Otherwise  results  would  be  of 
academic  interest  only.  The  method  presented  here  will  enable  the  designer 
to  quickly  select  a  system  with 'the  required  gain,  time  lag,  and  power  out- 
put from  the  servo  drive  and  install  it  near  the  optimum  position  relative 
to  the  airplane  reference  axis.  Flir^t  testing  tlJ^  system  from  that  point 
would  probably  yield  the  quickest  dependable  results  unless  an  smalog  com- 
puter were  available. 


i 
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■  CALCULATIONS 

The  lateral  equations  of  motion  are  given  in  reference  9.      In  all 

cases     Cy     ,     Cy   .   and    ^     were  assumed  tb  be  zero.     The  equations   including 

P  r 

the  product  of  inertia  were  used  only  in  case  I.     Otherwise  they  were  used 


as  follows: 


then» 


where 


(m  t)  -  C„   )q  -  C,;Z(  +  m  Df       =  0 

o          ^g  P         ^  ° 

p-                         0  ^ 

-C     6-  C  DCf  +   ( I  D-C  )D  4^         «  0 


b^A^  +  bjA^  +  bgA^  +  b^A  +  b      »A  =  0 


b .  '  m  I  I 
4         o  X  z 

b-  *  -(m  Cn     +  Cv  I  )I     -  m  I  C 

bo*CvC>,    I     +(mC(7     +C«.I)C       +mIC       -mC..C 

b,   =  -C^  C^  C       +  m  C^  C    .  -  m  C/)  C     -  C^C.    I  +  C^  C^  ,  C   . 
^  ^^^^%        ^^""^        °^0°P       ^^p^       ^p'^V^'^iZ^ 

b     »  C.(C^  C         -  C.    .    C      ) 
o         L'  -^p  n^         i?  y     np 

The  fastest  and  most  accurate  method  for  the  solution  of  the  stability 
quartic  was  to  use  a  combination  of  the  Lin  approximation  method  and  the 
Zimmerman  plot  as  explained  in  references  11  and  13.     With  the  Lin  method, 
values  could  be  obtained  that  would  lead  to  an  intersection  in  two  or  three 
points  with  the  Zimmerman  plot.      In  general,  there  were  two  real  roots  and  a 


\ 
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complex  pair,  although  where  a  large  value  of  +C  ^  was  used  there  were 

two  complex  pairs. 

Where  the  roots  were  A  , ,  A  „  and  a  +  i  P  » 

.695      , 

T,  A-,  -  seconds 

L/l  a 

P    =  —r~       seconds 

For  the  plots  in  figures  3  and  4,  where  each  point  represents  the 
solution  of  a  quart ic,  the  equations  for  ^a  %  o  ^   f\     were  used  where  the 
variables  were  C   and  C  .   determined  from  the  particular  value  of 
f  «  f  (a,P  )  for  that  point. 

Using  the  above  three  equations  of  motion,  a  general  solution  for 
forcing  function  inputs  of  A.  B,  and  C  respectively  was  made.  The  transient 
responses  took  the  following  forms: 

lf>  a^A^  +  a^A^  +  a^A  +  a^ 


where 


a,  *  Cm  I 
^  ox 


ao  »  AC     I    ■♦-  Bm  C  ,  -  C(m  C^     +  C^  I  ) 
2  np  X  o  n^         ^  o  ^^        Ip  x 

a     '  A(Cg  C       -  G     C.J-BC^.C       +  CC.  C 
1  ^p^l^        "pV        '''^\  V^P 

\  -  ^^L^np  -  ^Cl^J?p 


a- A     +  apA     +  a^A  + 


a 
o 
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where 


^ 


AI  I 
X  z 


-A(I  C  .  +  I  C.  )  -  Cm  I. 


X  n 


fi 


0  X 


a_  - 


J(tT 


A(C^  C  .  -  C  C,   )  +  B(C.I 

-BC.C  .  +  CC.C.  . 
L  ny     L  5-vy 

a„A  +  a/A  +  a 


m  C  )+Cm  Cfl 


.:2 


where 


Eo  "^  Bm  I 
2     0  z 


a,  «  ACq  I  -  B(m  C  .  +  I  C^  )  +  Cm  C.  . 
1     ^Q  ^  o  P*p  ^  ^B      o  ix^ 

a  «  -A(C^  C    -  C  C.  .  )  +   B(C„  C  ,  +  m  C   ) 
o      ^B  ^tf/    ^  ^H^  ^B  "4^    °  "e 

In  the  solution  for  the  disturbance  in  p,  ^(t)  became  p  ,  and  A, 
B,  and  C  were  Cy  ,  C^   ,  and  C  . 

In  the  response  to  the  aileron  step  function,  5"^)  became  o        and 

A,  B,  and  C  were  0,  C/i   »  and  C    • 

a         a 

A,  and  hence  the  roots  ^  i  ♦    p  *  ^'^^  a  T  ip  are  known. 

The  solution  for  the  typical  lateral  transient  response  takes  the 


form 


^iV  .  M 


g(t)=C^  +  C^i  ^  ■*-  ^2^    +  Kje.*^  sin(pt  +  F) 
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By  use  of  the  Heaviaide  Expansion  Theorem, 
C 


a 


a- A     +  a,,X     +  a,A     +  a 
o     -  or        2  ^r        1     r        o 

\f     ^  _  _ 

^         A   (4b.X     +  3b- A^  +  2b„A     +  b, ) 
r^     4    r  3     r  2    r        1 

Then  from  end  conditions 

-(C     +  C.  +  C   ) 
^  sin  F 

-P(C     +  C.   +  Cj 


tan  F 


o         1        Z' 


^5-^A->^2^2^^^^o^^l^^2) 


where  C„  represents  the  initial  condition  of  ^      ,  P  ,  or  0  at  t  «  0. 

^3 
Cj.  =  r-  .  An  initially  relaxed  system  is  assumed. 

4 

To  find  the  frequency  response  for  the  critical  time  lag  calculations, 
standard  servoraechanlsms  methods  were  used.   It  is  believed  that  the  best 
method  of  reducing  the  data  is  outlined  in  reference  8,  and  is  as  follows; 

For  the  airplane, 

3      2 

^)  ^°  ■•'  ^2°  +  a,  D  +  a 


^r    b^D^  +  bjD^  +  bgD^  +  b^D  +  b^ 

Using  the  general  equations  developed  previously  with  forcing  func- 
tions A,  B,  and  C  as     C^       ain  uJt,     Qn       sin  ^  ♦  an<^     C         sin  ujt     the  a's 

T  T  "t 

can  be  quickly  calculated. 

Then  K.  ,  the  steady-state  frequency  response  of  the  airplane,  can 
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be  expressed 


K 


A 


5r 


^1  ^  ^^1 


5, 


where  the  substitution  D  =  ia)  Is  made. 


-aoCO  4-  a 
d  o 


X„  »  u)^  -  b^oo  +  b 


-a»a)  +  a,a) 


'^' 


Yg  «  -bjo)  +  b^o) 


-1  ^1 
^,  «  tan    Y^ 


^r>  «  tan 


-1  h 


Then  for  each  frequency  the  magnitude  of  the  airplane  gain  or  ampli- 


tude is 


JL 
^r 


X^  cos  0^ 
Xg  cos  0^ 


and  the  phase  angle  is 


0  -  ^1  •  ^2 
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T/IBLE  I 

Principal  Dimensions 

S  «  250  sq.  ft. 

b  «  35.25  ft. 

c^  «  7.45  ft. 

1^  «  13.45  ft, 

3  .  »  36  sq.ft. 

S  =  8.54  sq.  ft. 

c  .   =  6.75  ft. 
vt 

b  .   »  7.42  sq.  ft. 
W    «  12,000  pounds  plus  fuel 
Total  control  throw  for  the  rudder,  25+1  each  way^ 
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TABLE      II 

Condition  I                    Condition  II 

Condition  III 

CR  max     clean 

50^  NRP     clean 

Landing 

h 

35,000 

10,000 

sea   level 

V 

660 

259 

183 

M 

.679 

.323 

.134 

P 

.736(10)"' 

1.756(10)'' 

2.378(10)'^ 

W 

13.500 

13,500 

12,500 

M- 

64.6 

27.1 

20.0 

^ 

160.5 

59.2 

40  ;o 

fflo 

6.8 

7.38 

7.71 

b/27 

.0267 

.0681 

.0964 

Cl 

.337 

1.18 

1.20 

OC 

2.9 

8.7 

7.8 

Q 

3.63 

3.63 

4.07 

\ 

-  .73 

5.07 

3.73 

U    . 

30 

30 

30 

3 

-27.1                           -21.3 

-22.2 

^x.. 

7778                               7778 

1 

7822 

^Ho 

24244 

24244 

24921 

^X 

.0055 

.0149 

.0220 

l£ 

.0172 

.0465 

.0688 

^xi 

.00015 

-   .0028 

"  .0030 

K 

.5 

.817 

1.0 

Gp 

-  .546 

^.400 

^  .400 

0,x 

.132 

.174 

.231 

G, 

.             -   .119 

-  .100 

-  .064 
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Gi, 


'^€. 


Cj? 


A      ! 


0 


rt 

Gnx: 

% 
P 


L. 


1    ACn^ 


n 


H' 
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TABLE   II      (continued) 


.0103 

-  .0146 

0OO35 
+  .0380 

-  .230 
.135 

-  .0573 

-  .0062 
.0010 

-  .671 

2.26 
3.15 
1.26 
spirally  stable 


.0043 

.0023 
.082 

-  .0272 

-  .0386 

.0119 

.0223 

-  .0350 

-  .0520 

-  .175 

-  .125 

.128 

.135 

-  .0444 

-  .0410 
.015 

-  .0119 

-  .0121 

-  .0024 

-  . 0050 

-  .562 

-  .344 
.115 

3.56 

^.13 

7.93 

8.28 

1.94 

2.05 

31.5 

15.1 

With  Yaw  Damper  Operative 


-  .0255 
.0137 

-  .0315 
.0147 
2.44 

.323 
.13 


-.0341 

.0129 

-.0450 

.0115 
4.28 

1.  35 

.31 


-.0379 

.0155 

-.0500 

.0105 
4.73 
2.02 
.43 
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TABLE      III 


Sln£:le  Degree  of  Freedom 


Condition  I 

:  iz 

.0172 

Gra» 

- .0061 

Gn^ 

.135 

(jJo 

2.805 

CJ,i 

2.80 

^ 

.064 

P 

2.24 

T'/z 

3.88 

C/. 

1.72 

K  for 
Q-  .5 

.26 

Gy,  If 
K=  1 

.13 



Condition  II 


Condition  III 


.0465 
-  .0119 

.128 
1.661 
1.658 
.077 
3.78 
5.66 
1.43 

.51 
.304 


.0688 
-.0121 

.135 
1.404 
1.40 
.063 
4.48 
7.90 
1.76 

.76 

.402 
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